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Abstract 



When computing the properties of reactions involving unstable charged particles care 
has to be taken to use a gauge invariant amplitude. In this talk we present methods 
to (automatically) obtain such an amplitude, both at the tree level and in one-loop 
^ ' calculations, using only a minimal number of diagrams. The numerical difference 

^■f^ ^ with gauge variant methods commonly used will be discussed for the case of tree 

, level W pair production. 
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G^ ■ 1 Introduction 



^ • Many of the cross sections one would like to evaluate with the methods described 

I I elsewhere in these proceedings involve the production and decay of unstable parti- 

Oh. cles. The W and Z bosons both have a sizable width T (F/m « 1/40), and the same 

' may be true for the top quark and Higgs boson (depending on their masses). For 

these particles the well-known narrow width approximation will not be sufficient 
for many purposes. Including off-shell effects, however, involves a resummation of 
\ some of the higher order graphs, which means that care has to be taken to keep 

, the matrix element gauge invariant. This last point is not only of academic interest 

' as gauge breaking terms are often much larger numerically than the gauge invari- 

ant result. Therefore even suppression factors F/m may not be enough to avoid 
inaccuracies. 

The Dyson resummation of higher order graphs in the propagator, which leads to 
a finite width, is clearly only needed when the unstable particle can kincmatically be 
on its mass shell. Indeed, for an unstable particle with space-like momentum < 
the imaginary part of the self energy is zero, hence no finite width should be used. 
One can therefore not use a complex mass in the Lagrangian, even disregarding the 
problems with the standard model mass generation mechanism which would follow 
(for instance a complex value for sin^ Ow)- It is a purely kinematical problem which 
particles should be given a finite width and which not. 

The amplitude can thus be divided into a set of resonant diagrams, which contain 
a divergence within the allowed phase space as an unstable particle is taken on-shcll, 
and non-resonant graphs without such divergences. In principle the contribution 
of the non-resonant diagrams is suppressed with respect to the resonant ones by 
T/m cx a. We first discuss ways to only use the resonant tree graphs to get the 
dominant behaviour, next increasing the accuracy of the result by including the non- 
resonant tree graphs and the resonant one- loop graphs. A more detailed discussion 
will be given in a forthcoming publication The tree level results are based on 
Ref. il. 
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2 Resonant Tree Graphs 



In this section some schemes will be given to deal with the resonant tree graphs. 
For simplicity these are first discussed for a single unstable particle with momentum 
p and mass m; the extension to more particles will be given later later. We assume 
that the diagram generator has already sorted out which diagrams ones are resonant 
and which are not. Without resummation the full (gauge invariant) amplitude for 
this process therefore has the form 

M = 4^ + N{p'). (1) 

The essential variable here is the virtuality of the unstable particle. The other 
kinematical variables are assumed to be CMS angles, which are independent of p^ . 
Unfortunately, this amplitude gives an infinite cross section when the integration 
region includes the pole at p^ = vn? . 

The simplest approximation one can make is to completely factorize the process 
in production and decay: the narrow width approximation. This fails to take into 
account terms of OiV); the corresponding dimensionless parameter will usually 
be r/m (about 1/40 for the W), but can be much larger near the threshold for 
production of the unstable particles. 

The most important of these OiT) corrections are included by simply using the 
on-shell expression for the matrix element, but treating the kinematics and unstable 
particle propagator off-shell. This off-shell propagator is derived by resumming 
the one-loop corrections to the propagator, giving rise to a simple Breit-Wigner 
propagator l/(p^ — w? -I- imV) with T defined by the relation mV = Imn(m^) in 
terms of the self energy 11 evaluated at the real mass. We refer to this procedure 
as the "resonant" scheme. This procedure obviously respects gauge invariance and 
will serve as the base line to compare the others with. 

The next step is usually to also evaluate the matrix element off-shell, we refer to 
this procedure as the "off-shell" method. This also gives an answer below threshold. 
In the case of a charged resonance however, this procedure violates U{l)em gauge 
invariance, even if the width is given by the physical (on-shell) quantity, which is 
gauge invariant. The reason is that the original resonant graphs are not the only 
graphs which lead to the particular final state: non-resonant graphs have to be 
included. 



3 Non-resonant Tree Graphs 

One reason to include the non-resonant graphs is to conserve gauge invariance. 
Another is that they can be sizable due to some other near-divergence, such as 
large collinear logarithms. The "naive" way to include these non-resonant graphs 
is to just add them to the resonant terms with a complex mass. However, as the 
amplitude was gauge invariant for F = it follows that it must be gauge variant 
for F > (if the resonant graphs are not separately gauge invariant). 

An obviously gauge invariant procedure to include the non-resonant graphs was 
introduced by Zeppenfeld et al a. The price to save gauge invariance in this 
"overall" scheme is an incorrect treatment of the non-resonant contribution close 
to mass shell. However, it can be argued that the difference is of higher order in 
this region of phase space. We thus have a prescription which is correct to leading 
order both on resonance and away from it, but these contributions are formally of 
different order in a. 

Another gauge invariant way to include the off-shell effects is to systematically 
separate orders in F|^. As the residue at the pole p^ — is gauge invariant 
one can add the finite width in the first term without breaking gauge invariance. 
This corresponds to adding and resumming only the gauge invariant part of the 
propagator corrections. In this scheme the cross section is given as the sum of the 
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Table 1: The six different schemes to treat unstable particles at tree level. 



resonant cross section plus 0{r) corrections, which are both gauge invariant. It is 
thus a natural starting point for higher order corrections. 

However, this "polescheme" also has some undesirable properties. The first 
(resonant) term has a discontinuity when the threshold for the production of the 
unstable particle is crossed. Approaching from below, one even encounters the 
original (non-resummed) singularity in the propagator. The accuracy of this scheme 
is thus doubtful around threshold. 

The difference between the last three schemes is of order for total cross 
sections and singly differential distributions, but of order T for doubly and higher 
differential distributions. This follows from the occurrence of a Levi-Civita tensor 
in the 0{T) terms. A summary of all schemes is given in table 

The extensions to the case where there are more unstable particles are largely 
straightforward. There now appears a hierarchy of n-fold resonant diagrams. The 
first three methods just use the maximally resonant ones. In the "polescheme" one 
only needs these and the ones where only one particle is non- resonant; the "overall" 
scheme uses all diagrams. The threshold behaviour of the 'polescheme" now is 
worse than before. However, in this area the worth of any fixed-order calculation is 
doubtful because of possible bound state effects. 

The numerical influence on the total cross section is illustrated in Fig. |l|, where 
the total cross section for 77 £~DgW~^ is plotted, with the considered stable 
and a photon beam derived from backscattered electrons [|| . The physically relevant 
off-shell effects will be twice as large. To avoid the collinear divergence in the final 
state a minimum transverse momentum p±{£) > 0.02-^3 is demanded; this should 
also give an indication of detector acceptance. It can be seen that the narrow width 
approximation misses some important effects. The graph for the offshell scheme is 
misleading, in that we used the formula ef = —g'^'^ which is not valid when 
the photon couples to a non-conserved current. It does show that one has to be 
careful when breaking gauge invariance. One should also note that the differences 
between the various schemes are expected to be much larger in doubly differential 
distributions. 
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4 Resonant One-loop Corrections 

When computing the 0{a) corrections one is normaUy only interested in the reso- 
nant part. We wiU thus have to identify which diagrams contribute to this, and then 
to find a gauge invariant way to isolate just the resonant part. A generalization of 
the "polescheme" is proposed for this last step. In the case of a neutral resonance 
(like the Z boson) the isolation step is not necessary, as the resonant diagrams al- 
ready form a gauge invariant subset; the only bone of contention left in this case 
amounts to the proper definition of the mass of the unstable particle. 

There are three classes of resonant diagrams. The first one contains the correc- 
tions to the resonant propagator. A gauge invariant part of these are resummed to 
give the finite width; the rest (vanishing at = m?) has to be treated perturba- 
tively. The second group, the factorizable diagrams, comprises the corrections to 
either production or decay of an unstable particle; these diagrams obviously retain 
the original resonant propagators. Added to this, charged resonances will give rise 
to logarithms log(p^ — rn?) as the virtuality and finite width of the unstable particle 
regulate the infrared divergence. The third class of non-factorizable diagrams is 
given by the set of infrared divergent diagrams with the photon (or gluon) spanning 
one or more resonant propagators. Simple power counting arguments show that 
no other non-factorizable diagrams can give a resonant contribution. One would 
suspect that even the contributions from these diagrams may cancel against the 
corresponding soft Bremsstrahlung diagrams. As the hard Bremsstrahlung is eas- 
ily seen to be non-resonant this would imply the intuitively satisfying result that 
initial-final and final-final state interaction vanish in the narrow width approxima- 
tion. 

However, this is not the case. As an example consider the doubly resonant five 
point function shown in Fig. ^ Using the Schouten identity one can decompose it 
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Figure 2: A double resonant five point function occurring in W pair production 

in five four point functions as 

= 2A^ (a?s^o(7) + A^psDoif^) + A^"Do(W^-) + A^/ DoiW+) + A^,Do(e)) 

(2) 

where Dq (p) denotes the scalar four point function with the propagator correspond- 
ing to particle p missing and the A's denote determinants of external and internal 
momenta Q. The overall factor — I/2A5S is quadratic in the virtualities p'j. — vn? 
of the W bosons and thus gives rise to the double resonance. The -Do (7) and A^^ 
are both finite, whereas for the other combinations the determinant is linear in the 
virtuality and the four point function linearly divergent. 

The corresponding scalar Bremsstrahlung integral can be decomposed similarly, 
except that the equivalent of the -Do (7) never appears. The other four terms exactly 
match the corresponding terms in the virtual diagram, leaving this one doubly 
resonant term. There is no symmetry which would cause it to disappear, even 
when integrated over phase space, and there are no other graphs which will give 
rise to the same four point function. It therefore is a double resonant term of order 
a which is not suppressed by additional factors T/m. The failing cancellations can 
be traced back to the fact that, in the fimit F <C A (A is an infrared regulator, 
normally the photon mass) the linear and quadratic infrared divergences do not 
cancel. Only a numerical study can decide whether the non-factorizable diagrams 
give rise to measurable effects. 

The diagrams we have selected as resonant diagrams in general do not form a 
gauge invariant subset. One can select a gauge invariant part by taking the residue 
at the pole, as was done in the previous sections. Unfortunately, this residue is ill- 
defined due to the occurrence of the divergent logarithms. (These also do not cancel 
between the virtual and Bremsstrahlung graphs.) These logarithms only occur in 
the scalar functions however. As the coefficients of the scalar functions are gauge 
invariant, one can evaluate these any way one likes. 

The following recipe therefore serves as the extension of the "polescheme" for 
one-loop (and soft Bremsstrahlung) diagrams: 

• The reduction of the diagram to scalar functions should be done on-shell (if 
the unstable particle is charged). 

• The scalar functions can be evaluated in any convenient non-divergent way; 
on-shell or off-shell. This also allows one to take into account the finite width 
in non-divergent but non-analytical logarithms like the Coulomb singularity 
and threshold effects. 

This schemes have been implemented in the FF library 

5 Conclusion 

In this talk we have presented a scheme to systematically include off-shell effects 
in amplitudes including unstable particles, explicitly keeping gauge invariance. Un- 
fortunately it breaks down near the threshold for the production of the unstable 
particles. In this region, however, the occurrence of bound state effects makes any 
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fixed-order perturbation series suspect. The numerical difference of this procedure 
with the (gauge variant) ones commonly used are small in the tree level total cross 
section, at least in the case of W pair production and the unitary gauge. Larger 
differences are expected for doubly differential cross sections. 
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